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statistics of the random sample (X±, Yi), {X%, Y%),. . . , (-Xjv(t)i ^jv(t)) which is then used to obtain asymptotic 
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specific insurance settings. As a by-product we obtain a stochastic representation of a m-dimensional A- 
cxtrcmal variate in terms of iid unit exponential random variables. 
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1 Introduction 

Let {(A„, Y n ), n > 1} be bivariate random claim sizes with common distribution function F arising from an 
insurance portfolio. In a specific insurance context, X n may stay for instance for the total claim amount 
related to the nth accident, and Y n for the corresponding total expense amount. Let further {N(t),t > 0} 
be a stochastic process which counts the number of claims that occur in the time interval [0,i],i > 0. So 
we observe (X\, Yi), . . . , (Ajv(t), Yv(t)) claims up to time t > 0. Denote by Aj.jy(t), Y i:N ( t ), 1 < i < N(t) the 
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corresponding ith lowest claim in the above random sample taken component-wise. 

Especially in reinsurance applications, the largest claims -X"./v(t)-i+i:JV(t)j YN(t)-i+i-.N(t)i i = 1, ..-,771 are of 
particular importance. For instance, consider the following reinsurance contract introduced by Ammeter 
(1964) 

Sl(p,t) = -X"jv($ : jv(t) + ^N(t)-l:N(t) + 1" ^N(t)-p+l:N(t)j P>1 

which is known as the Largest Claim Reinsurance. Thus at the time point t a reinsurer covers the total loss 
amount S\ (p, t) which is the simplest linear transformation of the upper p largest claims observed up to time 
t. Another simple reinsurance treaty 

Sl(p,t) = ^j\r(t):iV(t) + ^N(t)-l:N(t) + 1" ^N(t)-p+2:N(t) ~ (P ~ ^)^N{t)-p+l:N(t)> P > ^ 

is introduced by Thepaut (1950). This treaty is known in actuarial literature as the ECOMOR reinsurance 
treaty (for further details see e.g. Teugels (1985), Beirlant et al. (1996), Embrcchts ct al. (1997), Rolski et al. 
(1999), and Mikosch (2004)). 

In general one can define the total loss amount by 

p 

S\ip,t) = y^g*(^jy(t)-j+i:iv(t)), 
j'=i 

with g* , 1 < j < p real measurable functions. Typically, the prime interest of a reinsurer is the calculation 
of the pure premium E{S\(p, t)}. This can be clearly expressed as $^j=i E{g*(X N ^_j +1:N ^)} supposing 
additionally that fl^CXjvft)— j+i:;v(t))) 1 < j < P have finite expectations. It is often assumed in actuarial 
applications that N(t) is independent of claim sizes. With that assumption we have for each j = l,...,p 

oo 

E{g*(X N{t) _ j+1:N(t) )} = E{E{g*(X N(t) _ j+1 .. N(t) )\N(t)}} =J2E{g*(X k _ j+1:k )}P{N(t) = k}. 

k=l 

Essentially, in order to compute the pure premium, we need to calculate the expectations of the upper order 
statistics. Further we need to know the distribution function of the counting random variable N(t). Next, 
since we consider a bivariate setup, let us suppose further that a similar treaty 

S2(q,t) = ^2gj*(Y N ( t )-j+i:N(t)), ?>1) 
J=l 

with g** , 1 < j < q real measurable functions, covers the risks modelled by Y\, Y2, ■ ■ ■ , Viv(t)- 

In order to price both treaties mentioned above, the reinsurer needs to have some indications concerning the 
distribution function of the total loss amount (Si(p, t), ^(g, t)). For instance, if in particular the standard 
deviation (variance) premium principle is used, then an estimate of the standard deviation is required. Further 
the dependence between Si(p, t) and S2{q,t) needs to be quantified. 

An asymptotic model for both treaties can be regarded as a good candidate to overcome the difficulties 
in specification of the model (distribution assumptions for the claim sizes or assumptions on the first and 
second moment of X\,Yi). The idea is to let t — ► 00 and to investigate the joint asymptotic behaviour of 
(S 1 (p,t),S 2 (q,t)). 

Relying on extreme value theory Embrechts et al. (1997) (see Example 8.7.7 therein) derives the asymptotic 
distribution of both the Largest Claim Reinsurance and the ECOMOR Reinsurance treaties when N(t)/t — > 
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A G (0,oo) in probability. In Hashorva (2004) the asymptotic limiting distribution for (Si(p, t), S 2 (q, t)) a 
bivariate ECOMOR reinsurance treaty is obtained when the marginal distributions of F are in the max-domain 
of attraction of the Gumbcl distribution, imposing further the iid (independent and identically distributed) 
assumption on the claim sizes. 

With main impetus from the afore-mentioned results, we consider in this article special reinsurance treaties 
with g*j,g*j* simple linear functions and claim sizes which can be dependent (thus dropping the iid assumption). 
In Section 2 we first deal with the joint distribution function of randomly indexed upper order statistics; 
application to reinsurance and details for the iid case are presented in the Section 3. Proofs of the results are 
given in Section 4. 



2 Joint limiting distribution of randomly indexed 
upper order statistics 

It is well-known (see e.g. Reiss (1989)) that the joint convergence in distribution under the assumption of iid 
claim sizes 

X n .. n -bi(n) Y n:n — b 2 (ri)\ /X„_ m+1:n — &i(n) Y n - m+1:n — 62(71) 



«i(n) a 2 {n) J \ Oi(n) a 2 {n) 

^({Xi,y 1 ),...,{X m ,y m )), ji^co (2.1) 

with {Xii 3^i), i < m random vectors holds for all integers m € Af and given real functions flj(t) > 0, bi(t), i = 
1, 2 iff the underlying distribution function F of the random vector (X\, Y\) is in the max-domain of attraction 
of H, the max-stable bivariate distribution function of (Xi,yi), i-e. 



lim sup 



F t (a 1 (t)x + b 1 (t),a 2 {t)y + b 2 (t))-H(x,y) = (2.2) 



holds. See below l|2.6|l for the joint distribution of X\, . . . , X^, k > 1. For short we write the above fact as 
F G MDA(H). 

We note in passing that the standard notation —>■, "—^ which we use throughout in this paper mean con- 
vergence in distribution, convergence in probability and almost sure convergence, respectively. 

Actually l|2.2|l implies that for the marginal distributions of F (denoted here by Fi,F 2 ) we have Fi € 
MDA(Hi), i = 1,2, where the standard extreme value distribution function Hi is either of the following 

$ai(x) = exp(-x~ ai ), Qfj > 0, x > 0, 

or 

* Ql (x) = exp(-|x| Ql ), cti > 0, x < 0, 

or 



A(x) = exp(— exp(— x)), xGTZ, 
i.e. unit Frcchct, Wcibull or Gumbel distribution, respectively 
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For iid claim sizes, if N(t) A oo as t — ► oo and N(t) is independent of the claim sizes for all t large, then by 
Lemma 2.5.6 of Embrcchts ct al. (1997) the asymptotic relation l(2.2|l implies for any i > 1 

( X N{t) _ i+1:N ( t) - h(N(t)) Y N(t) _ i+1:N{t) -b 2 (N(t))\ d 

{ ^um ' ) -* {XM t ~* 00 - 

A different situation arises when transforming with aj(t), bi(t) instead of the random functions 
a,i(N(t)),bi(N(i)),i = 1,2. So, if we assume further that 

M A Z, t^oo (2.3) 

holds with Z such that P{Z > 0} = 1, then it follows (along the lines of Theorem 4.3.4 of Embrcchts et al. 
(1997)) that for all i > 1 

'XN{t)-i+l:N(t)—h(t) YN(t)-i+l:N(t) - h(t)\ d 

^ ' ) ~* {XiMl 

holds with (X*,y*) a new bivariate random vector. 

For the univariate case Theorem 4.3.4 of Embrechts et al. (1997) shows the distribution functions of X* . 
Proposition 2.2 of Hashorva (2003) gives an explicit expression for the joint distribution of (X^, . . . , X*). Our 
next result is more general. We consider the bivariate setup allowing claim sizes to be dependent, and further 
instead of 1(2.3(1 wc assume the convergence in distribution 

Z, t->oo, (2.4) 
with Z almost surely positive. Clearly the above condition is weaker than l|2.3|l . 

Proposition 2.1. Let {(X n ,Y n ),n > 1} be bivariate claim sizes with common distribution function F inde- 
pendent of the counting process {N(t),t > 0} for all t > 0. If condition (12.4(1 is fulfilled with Z positive and 
non-zero and further (|2.1fl holds for some fixed m <E Af with constants a,(t) > 0, bi(t),i = 1,2 that satisfy 
(|2.2|) . then we have for t — > oo 

^N(t):N(t) — h{t) yiv(t):JV(t) ~ b 2 (t) \ f ^-N{t)-m+\:N{t) ~ bl(t) Y N ( t )_ m+1:N ( t } — b 2 (t) 



ai(t) a 2 (t) / ' ' \ a 2 (t) 

S ((z^x 1 + 6 1 ,z^y 1 + 5 2 ),...,(z^x m + 6 1 ,z^-y m + 5 2 )y (2.5) 

with 7, := l/oti,-l/oti,0 if F, & MDA($ ai ), MDA($ ai ), MDA(K), respectively. Further 5 i := InZ if 
Fi £ MDA(A) and otherwise for i=l,2. 

In the above proposition we do not assume explicitly the independence of the claim sizes. For iid claim sizes 
(recall that (|2.2() is equivalent to we obtain immediately: 

Corollary 2.2. Let {(X n ,Y n ),n > 1} be iid bivariate claim sizes with distribution function F independent 
of N(t),t > 0. If condition l|2.2|l is satisfied and further (|2.4|) holds with Z almost surely positive, then (12.51) 
holds for any m > 1 with 5i,ji,i — 1,2 as in Provosition \2. 1\ 

Remarks: (i) If {X n ,n > 1} are iid with distribution function F 6 MDA(H), with H a univariate extreme 
value distribution, then the joint density function of (X±, . . . ,X m ),m > 1 is given by (see e.g. Embrechts et 
al. (1997) p. 201) 

m H'( ■) m 
h m (x) = H(x m ) Y\ -777— y, with xi > x 2 > ■ ■ ■ > x m , J H ( Xi ) 6(0,1). (2.6) 

i=i n[ - x% > i=i 
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Referring to Embrechts et al. (1997) the random vector (Xi, . . . ,X m ) is called a m-dimensional if -extremal 
variate. If if is the unit Gumbel distribution, then the following stochastic representation (see Theorem 7.1 
of Pakes and Stcutcl (1997)) 

(Xt - X i+ i) = (—) , k>l (2.7) 

V / i=i fe V i / i=\,...,k 

holds with Ei,i > 1 iid unit exponential random variables. The standard notation = means equality of 
distribution functions. See ((3.12(1 for a stochastic representation of (X\, . . . ,X m ). 

(ii) In the above proposition Z is independent of (A^,[Vi), 1 < i < m. This follows immediately recalling that 
the claim sizes are independent of the counting process N(t),t > 0. 

For the more general class of strictly stationary random sequences conditions for joint weak convergence 
of upper order statistics are available in literature. For a-mixing stationary sequences (univariate case) 
the mentioned weak convergence is discussed in Hsing (1988). In Theorem 4.1 and Theorem 4.2 therein 
necessarily and sufficient conditions for joint weak convergence are derived. Moreover, it is shown that the 
limit distribution should have a specific form. Novak (2002) gives a simplified version of the limit distribution. 
A variant of Theorem 4.2 of Hsing (1988) can be found in Leadbetter (1995). Convergence results for dependent 
bivariate random sequences arc obtained in Huslcr (1990,1993). 



3 Joint Asymptotic Distribution of S\(p,t) and S2 

In this section we investigate the asymptotic behaviour (t — > oo) of (S±(p, t), S 2 (q, t)) defined in the intro- 
duction. In general some restrictions on the choice of the functions g*, I < j < p,gj*, 1 < j < q should be 
imposed. Obviously the total expected loss should be non-negative. Tractable (simple) functions which we 
consider here are 

g*(x) = kjtx, l<j<p, and g**(x) = k j2 x, 1 < j < q, (3.8) 

with kjx,l < j < p,kj2,l < j < q real constants. Put throughout in the following c\ := X^=i^ji> c 2 : — 
y,|_, kj2- Based on the previous results, we consider now the joint asymptotic behaviour of the reinsurance 
treaties. 

Proposition 3.1. Let {(X n ,Y n ),n > 1}, {N(t),t > 0},Z and ji,5i,i = 1,2 be as in Proposition^^ and 
for p,q G let g*, 1 < j < p, g** , 1 < j < q be as in 1)3.8(1 . If condition l|2.3(l is fulfilled and 1(2.1(1 holds for 
some m > max(p, q), then 

( 5i(p,t) - fri(f)ci S 2 {q,t) - b 2 (t)c 2 \ d / 71 ^ x 7 , ^ x \ , , 

{ aW) ' aW) ) [ 211 ^^+ c ^ 272 ^ j2 y j + c 2 6 2 j (3.9) 

holds as t — > oo. 

Remarks: (i) In the above proposition there is no restriction imposed on the constants kji. Neither need 
the claim sizes be positive. Typically, in reinsurance claim sizes are assumed to be positive, and further, 
some restrictions have to be imposed on the constants kji. Implicitly we may require that these constants 
are such that the expected loss (at any time point) for both Si(p,t), S 2 (q,t) is non-negative. A more explicit 
assumption would be to suppose that both c\,c 2 are positive. 
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(ii) Proposition 12.11 and Proposition 13.11 can be shown with similar arguments for a general multivariate 
setup, i.e. for claim sizes being random vectors in lZ k , k > 3. 

We discuss next the bivariate counterpart of Example 8.7.7 of Embrcchts et al. (1997), which was our starting 
point. We correct a missing constant in an asymptotic result therein. 

In order to keep things simple, we impose throughout in the following the assumptions of Corollary 12.21 
considering therefore iid claims sizes independent of the counting process N(t). 

Note in passing that in Example 8.7.7 of Embrechts et al. (1997) the claim sizes (univariate setup) are iid 
being further independent of the counting process, which satisfies (|2.3|) with 2 = A £ (0, oo) almost surely. 

To this end, suppose for simplicity that both marginal distributions F\, F 2 of F are identical and let F-[~ , F 2 ~ 
be the generalised inverse of F\ and F 2 , respectively. 

Case a): If F\ e MDA($ a ),a > 0, then we may take b\{t) = b 2 (t) = and 

0l (t) = a 2 (t) := M{x 6 K : F x (x) > 1 - 1/t} = Ff(l - 1/t), t > 1 

(see e.g. Resnick (1987), Reiss (1989) or Embrechts et al. (1997)). So we obtain for p,q e Af and kji real 
constants 

as t — > oo. Embrechts et al. (1997) proves for the ECOMOR case 

Si(p,n) d x i/o. <v v \ 

— ► A ' y i(Xi-X i+1 ), n^oo 

oi (n) ^ 

which follows immediately from i|3.10[l putting Z = A > 0. 

Case b): When F\ £ MDA(fy a ),a > holds, then we may take bi(t) = b 2 (t) := > 1, with cj := sup{a; : 
Fi(x) < 1} the upper endpoint of the distribution function F\ which is necessarily finite and define 

ai(i) = a 2 {t) := uj -Ff(l - 1/t), i > 1. 

Thus we have 

/ Si(p,t) - ciu S 2 (q,t) - c 2 uj\ d 



ai(t) a 2 (t) 



J'=l 3=1 / 



Case c): If Fi e MDj4(A), then we put for i large 

bi(t) = b 2 (t) :=Ff(l-l/t) 

and 

oi(t) = oa (t) := / [1 - Fi(s)]/[1 - Fi(6i(t))] da. 

Jbl(t) 

Thus we have that the right hand side of i|3.9[) is given by 

/ p q \ 

kjiXj + Ci In Z, ^ ^2^- + C2 In Z . 
\i=i j=i / 



G 



Next, we consider 3 examples. 

Example 1. (Generalised ECOMOR reinsurance treaty). Assume that the constants fey are such that 
Ci = C2 = 0. This is fulfilled in the special case of the ECOMOR treaty (see introduction above). 

Thus under the assumptions of Proposition ^. II we have for Fx, F% G MDA(A) the convergence in distribution 
/ S 2 (q,t) \ d (s^ u y, ST U v\ , _ 

For the ECOMOR treaty Embrcchts et al. (1997) obtains with Z = A almost surely 



p-i 



Si(p,n) d^ 
oi (n) 



oo, 



with Ej, 1 < j < p — 1 iid unit exponential random variables, which follows also by the bivariatc result above 
(recall l(277)l1. 

There is a remarkable fact in the above asymptotic result, namely the random variable Z does not appear 
in the right hand side of the asymptotic expression. This is not the case in general for Fi,i = 1,2 in the 
max-domain of attraction of Frechet or Wcibull. 

Example 2. (Asymptotic independence of the components of the maximum claim sizes.) As mentioned in the 
introduction, for the asymptotic considerations, we are not directly interested in the joint distribution function 
F of the claim sizes, but on the limiting distribution H. In some applications, even if F is not a product 
distribution, it may happen that H is a product distribution, meaning H(x,y) = Hi(x)H2(y),Wx,y G 7Z. 
This implies that Xi is independent of 3^ for any i,j > 1 and thus [Si(p, t) — bi(t)c\]/ a\{t) is asymptotically 
independent of [S2(q,t) — &2(i) c 2]/ a 2(£) as t — > oo. So the asymptotic distribution of each treaty can be 
easily calculated using l|2.6|l . We note in passing that there are several known conditions for asymptotic 
independence, see e.g. Galambos (1987), Resnick (1987), Reiss (1989), Falk et al. (1994), Hiisler (1994). 

Example 3. (F±, F^ with exponential tails and N(t) Poisson). Consider iid claim sizes with joint distribution 
function F which has marginal distributions tail equivalent to the unit exponential distribution function, i.e. 
lirrix-too exp(x)[l — Fi(x)] = l,i = 1,2. It follows that F\,F% G MDA(A) with constants a\(t) = a^{t) = 
1, foi(t) = b2(t) = In t,i > 0. As in the above example, assume further that the distribution function H is a 
product distribution and N(t) is a (homogeneous) Poisson process with parameter A > independent of the 
claim sizes. We have thus N(t)/t < ^»' A and for any p, q G Af 

^Si(p, t) - ci hit, S 2 (q,t) - c 2 Intj -4 kjiXj + ci In A, kj%yj + c% In , t — > oo, 

with Xj independent of 3^- Recall c\ := Y^j=i kji> c 2 := Ylj=i kj2- 

Borrowing the idea of Example 8.7.7 of Embrechts et al. (1997) we find an explicit formula for the right 
hand side above. Let therefore Eji,j > l,i = 1,2 be iid unit exponential random variables and put ku := 
Ylj=i kji/l, I > l,i = 1, 2. It is well-known that (see e.g. Reiss (1989)) 

n p 

E n - j+lm!i ) =(E^ ■ »>l,i = 1.2. 

/ j=l,...,n V, — 4 
1=3 



It is well-known that for large n 

Jt 1 

Inn = K + o(l) 



" 1 

1=1 
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where K is the Euler-Mascheroni constant. So we may write for n > p > 1 

p p n 

" / 



E fcji (Sn-j+im.i - Inn) = £ %i £ — ci In n 

i=l j=l i=j 

/ 



f 1 TP 1 

E fc u-% + c i £ ~ 1 i " Ci (^ _ Et)' n °°- 



On the other hand 



I 

i=i i=p+i i=i 



(E n -.n,i -Itlti,.. .,E n -j + i:n t i —]nn) — > (Xl,...,Xj), n — > 00, 



hence the continuous mapping theorem (see e.g. Kallenbcrg (1997)) implies 

p p 

— Inn) — > £ kjiXj, n — > oo. 
i=i i=i 

Thus we have the stochastic representation 

p p 



£%i^ = ^fca^i + ci £ ^L-^ + ciXp, (3.11) 



i=i 2=1 Z=p+1 



with Ki := K — Y?i=i 7)* — 1- We obtain thus proceeding similarly for the second treaty and recalling the 



i 

asymptotic independence assumption 



Si (p,t) — Cilnt, S 2 (q, t) — c 2 In t 

( ^ < ~ K ^ TP 1 y TP 1 

'£fcji% + ci[ £ 3-1 ,~ + if p + In a] , £ fej-a-Eja + [ E j2 .~ + Jf g + lnA 



VJ = 1 J=P+1 J 3=1 J=9+l ^ 

Remarks, (i) For a suitable choice of constants (|3.11() implies for any m > 2 

(x u ..., x m ) = (e u + £ + Ku ■ ■ ■ ' Smi + £ ^r^ +Km )> ( 3 - 12 ) 

1=2 l=m+l 

hence in particular we have for any i G M 

00 i 2 * 1 

E{X i } = l + K i , Var{^} = l+ E 72=^ + 1 -E72 < 



P 6 1 " P 



(ii) In Example 8.7.7 of Embrechts et al. (1997) the term — fclnfc in page 519 should be k(K - Y/ 3 =i })• 
A general result is given in the following proposition. 

Proposition 3.2. Let {{X n , Y n ), n > 1} be iid bivariate claim sizes with distribution function F. Assume that 
(|2.2|) holds with given functions a.i(t) > 0,bi(t),i = 1,2 and H{x,y) = exp(— exp(— x) — exp(— y)), x, y G 7Z. 
Suppose further that {N(t),t > 0} is a counting process independent of the claim sizes satisfying <|2.4[) with 
Z positive and non-zero. Then we have the convergence in distribution 

( Sijp, t)~bi{t)ci S 2 (q,t)-b 2 (t)c 2 \ 
I Oi (t) ' a 2 (t) j 
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(V oo „ q oo 

£ k^Ej! + In Z + K p ] , £ fc i2 S i2 + c 2 [ E + In Z + 

3=1 J=P+1 J i=l j=q+l 3 

(3.13) 

as t — > oo, iot£/i Ci, kji, K p , K q as defined above and Eji,j > 1,£ = 1,2 £id itm£ exponential random variables. 

4 Proofs 

Proof of Proposition I2TH Define in the following for t positive 

rrt /AT/ f « ._ (T fNfA S rp lN ( f \\\ -I X N{t)- 3 + l:N(t)-bl{N{t)) Y N{t) - j + 1:N{t) -b 2 (N(t)) \ 

T 3 (N(t)) . (T 3l (N(t)),T 32 (N(t))) , J 

and 

Z t := — — , a,i{t,z):= — -— , bi{t,z):= —r , z> 0,2 = 1,2. 

i ffli(t) <w) 

Let 0t,t > be arbitrary positive constants such that lim^oo z t — z G (0, oo). Assumption (|2.2|) implies that 

the positive norming functions a\,a 2 are regularly varying (see e.g. Resnick (1987)), hence using Theorem 

A3. 2 of Embrechts et al. (1997) we obtain 

lim di(t,Zt) = z 7i , i = 1,2, 

t — >oo 

with 7l := I/an, -l/a 4) if Fi G MDA($ ai ), F t G MDA{^ ai ) or F, G MDA(A), respectively. 
If F G MDA(^ a ) or F G MDi($ a ) with a > 0, then for large t we have b,(ts) - b % (t) = 0, i = 1, 2, Vs > 0. 
It is well-known (see Proposition 0.10, Proposition 1.1, Corollary 1.7 and Exercise 0.4.3.2 of Resnick (1987)) 
that if F G MDA(K) then lim t _, 00 (bi (ts) - bi(t))/a,i(t) = lns,s > holds locally unifromly on (0, oo). 
Consequently we have 

lim bi(t, Zt) = Sihxz, £ = 1,2, 

t — >oc 

with 5i = 1 if Fi G MFM(A) and S t = 0, otherwise. 

In view of (|2.4|) applying now Theorem 3.27 of Kallenberg (1997) we obtain the convergence in distributions 



ctj(t, Zt) — > 2 7i , and &i(i, Zt) — > <Jj lniJ, i— >oo,i = l,2, 



hence as i — > oo 



U t ,N(t) ~ (ai(t,Z t )Mt,Zt)Mi>Zt)Mt,Zt)) ^ {Zi\Z^,8i\nZ,8 2 \nZ) =: U z . 

Let KJj G 1Z 2 , 1 < j < m and m G (0, oo) 2 x 1Z 2 be given constants. Since the claim sizes are independent of 
N(t), t > we obtain by conditioning 

P{U t ,N(t) <u,T x (N{t)) <w 1 ,...,T m (N(t)) <w rn \N(t)=n} = l(U t>n < u)h(n), n>m,t>0, 

with h(n) := P{T\(n) < wi, . . . , T m (n) < w m }. Here a < b, a, b G 7Z k , k > 2 is understood component-wise 
and l(-) is the indicator function. By the assumptions 

lim h(n) = P{(X 1 ,y l )<wi,...,(X m ,y m )<w m }. 

n — >oo 

Using (|2.4|l we obtain N(t) A oo as t — > oo implying 

h(N(t)) ^ P{(x x ,y 1 )<w 1 ,...,(M m ,y m )< Wm }, t^oo. 
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If further u is such that P{U z = u} = then 

l(U tt N(t) <«)->• l(Uz < u), i^oo, 
consequently by Slutsky lemma (see e.g. Kallenberg (1997)) 

i(u ttm )h(N(t)) 4. i(u z <u)p{(x 1 ,y 1 )<w 1 ,...,(x m ,y m )<w m }, t^oo. 

Since l(U tt N(t)) a (N(t)),t > is positive and bounded (consequently uniformly integrable) we have taking 
the expectation with respect to N(t) and passing to limit 

(ut^TtiNit)), . . . , r m (2v(*))) 4 (u z , {Xx, yx), (x m , y m )) , *-><». 

Now, for any j < m, t > we may write 

'-XjV(t)-i+l:JV(t) - h(t) yW(t)-j + l:7V(t) — &2(*) 



ai(t) a 2 {t) 

oi(t, Z t )2ji(JV(«)) + 6i(t, Z t ), 5a Z t )T j2 (N(t)) + b 2 (t, Z t )^ 

hence the proof follows by the continuous mapping theorem. □ 

Proof of Proposition 13711 By the assumptions and Proposition 12.11 we have that (|2.5|) holds, hence we 
may write for t > 0,p, q € M using the continuous mapping theorem 

/ Si(p,t) -bi(t)ci S 2 {g,t) -b 2 (t)c 2 \ 

\ oi (t) ' as (t) y 

= (^^it^JVW-j+irivCi) -h{t)]/ai{t),Y, k 3A Y N(t)-j+i:N(t) ~ M*)]/ ^*) J 

\3=1 3=1 / 

/ P 1 \ 

V 3=1 3=1 / 

with 6i, 7j, i = 1, 2 as in Proposition 12. II thus the proof is complete. □ 

Proof of Proposition 13.21 The proof follows immediately using Proposition 13.11 and the result of Example 
3 recalling further the expression for the joint density function given in l|2.6[) . □ 

Acknowledgment: I would like to thank two Referees for several comments and in particular for correcting 
the result of Example 3 pointing out a missing term related to Euler's constant. 
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